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Abstract 

We analyze the impact of using 6-tagged samples in studying non-Abelian effects 
due to QCD in e + e~ — > 4jet events at -/s = M z o, using angular variable analyses 
and comparisons with e + e~ — > 3jet7 events. We find that QCD effects are largely 
enhanced in 6-quark samples with respect to 'unflavoured' ones, where energy- 
ordering is used to distinguish between gluon and quark jets. We show that the 
6-quark mass influences the angular distributions significantly and should not be 
neglected. 
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1 Introduction 



In recent years the experiments ALEPH, DELPHI, L3 and OPAL at LEP I have performed a 
number of measurements of the process e + e~ — > Z° — > hadrons in order to point out effects 
due to QCD [|TJ. Important results have been achieved. The strong coupling constant a s has 
been determined from jet rates and from shape variables || and both the flavour independence 
H and the running with y/s have been verified. Three- [|5] and four-jet M distributions have 
been studied and their behaviour agrees with QCD predictions calculated to second order in ot s . 
The colour factors, which determine the gauge group responsible for strong interactions, have 
been measured @. Abelian models alternative to QCD have been ruled out and the coupling 
of the QCD triple gluon vertex has been verified to be in agreement with QCD predictions ||. 
Concerning the latter tests, several variables sensitive to differences between QCD and 



Abelian models have been proposed P, [10], 111], O, Oj. One of the main differences is the 



predicted relative contribution of e + e~ — > qqq'q' events in e + e~ — > 4jet samples: about 5% in 
QCD but 30% in Abelian models 0]. The tests are preferentially based on angular correlations 
between jets, which are usually ordered in energy and where the two most energetic jets are 
"identified" with the primary quarks (i.e., from the Z° decay). Among these angular variables, 
some of the most widely used in experimental analyses are the modified Nachtmann-Reiter 
angle 6* NR [[[(], |TT| , the Bengtsson-Zerwas angle xbz [0, the Korner-Schierholz-Willrodt angle 



®ksw 10 an d the angle between the two least energetic jets 6 U . The distributions in these 
variables are quite different for e + e~ — > qqgg and e + e~ — > qqq'q' events. 

Another way of searching for evidence of effects due to QCD, and only partially exploited 
so far, is to use photon samples. In particular, in order to isolate the triple gluon vertex 
contribution, one can compare 4jet- and 3jet7-samples. This approach is based on the similarity 



of photon and gluon bremsstrahlung off quarks (T^j and has already been adopted in ref. [|T 



However, there is also an obvious difference between photons and gluons, which is that only 
photons can also be radiated off the initial state electrons and positrons. Therefore some 
care is needed to eliminate the distortion due to this Initial State Radiation (ISR) from the 
3jet7-sample before a direct comparison with the 4jet sample can be made. 

It is the purpose of this paper to study to what extent the techniques of flavour identifi- 
cation, that are rapidly being developed by the various experimental collaborations at LEP I 
|PT|| , turn out to be useful in recognizing effects of QCD from the analysis of fe-tagged 4jet- and 
3jet7-samples, following both the approaches described above. The most widely used methods 
to recognize 6-quark jets are probably the following: 

• reconstruction of semileptonic 6-decays by observing a high lepton; 

• lifetime tagging by detecting a secondary vertex; 

• reconstruction using kinematical "event shape" or "jet shape" variables. 

Their main features are summarized in refs. |T3L |T3. Recently, in addition to these and other 



conventional methods |20j also the identification of gluon and quark jets by means of neural 



networks has been proposed [£l]] 



In our opinion, there is an important motivation for this study. As many authors have 
pointed out, the angular variables described here are most useful for emphasizing the non- 
Abelian features of QCD if one distinguishes between quark and gluon jets and assigns the 
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four-momenta of the final states to the corresponding particles. If that is not possible, the 
best one can do is to order the jets in energyf]. However, with 6-tagging, it is at least possible 
to distinguish some quark jets, namely those originating from 6-quarks, from gluon jets. This 
opens the prospect of observing the non-Abelian structure of QCD much more clearly by 
selecting events containing 6-jets. Since this means that all final states contain at least two 
6-quarks, the effect of the the 6-quark mass becomes much more important than it was before, 
and it must be properly taken into account in the analysis. The required matrix elements have 
only recently become available |22, 23, p4|P|. 



The plan of the paper is as follows. In section |] we give details of the calculation, of the 
algorithms used in the phenomenological analysis and the numerical values adopted for the 
various parameters. Section |] is devoted to a discussion of the differential distributions in four 
angular variables, and section (|to the possibility of using 3jet7 samples. In section || we study 
the sensitivity of our results to the 6-quark mass, and we draw our conclusions in section |]. 



2 Calculation 

The Feynman diagrams describing at tree-level the reactions 

e + + eT -> q 1 + q 2 + q' 3 + q' 4 , (1) 

e + + e~ -> q x + q 2 + 93 + 94, (2) 

e + + e~ -> qi + q 2 + g 3 + 74, (3) 

where q^ = u, d, s, c and 6, are shown in figs. 1-2. In the present analysis we have computed the 
corresponding matrix elements with the same FORTRAN code already used in refs. [^, p4j] , 
which takes all masses and both the 7 and Z° intermediate contributions into account exactly. 
For all the details of the computation, as well as for the explicit helicity amplitude formulae, 



we refer to 23 



The matrix element squared \A4\ 2 for the four quark process ([!]) can be written as 

\M\ 2 = Cf q ' q '\M + \ 2 + C 9 J q ' q '\M-\ 2 , (4) 

with 

CT V = \n c c f (t f t (c f - \c A )), (5) 

M ± = Mi + M 2 + M 3 + M 4 ± S qq '(M 5 + M 6 + M 7 + M 8 ), (6) 

where Mi is the amplitude corresponding to the i-th diagram in fig. 1, Nc = 3 the number of 
colours, Cp = (Nq — 1)/(2N C ) = 4/3 and Ca = N c the Casimir operators of the fundamental 
and adjoint representations of the gauge group SU(N C ), and T F = 1/2 the normalisation of 
the generators of the fundamental representation. 

2 In Z° — * qqgg events the lowest energy parton is a gluon in 1=3 84% of cases, whereas the percentage of 
events in which the two lowest energy partons are both gluons is only « 53% |w| . 

3 An updated review on matrix element computations for multi-jet production in e + e~ reactions presented 
in the literature can be found, e.g., in the introduction of ref. |23|. 
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The matrix element squared for the two quark and two gluon process (0) , can be split into 



two gauge invariant parts as follows pq| : 

\M\ 2 = Cf"\M a \ 2 + Cf"\M b \ 2 , (7) 

where 

M a =J2 M h (8) 

i=l,6 

M b = Mi + M 3 + M 5 - M 2 - Mi - M 6 - 2 i[M 7 + M%], (9) 
Mi, i = 1, ...8, corresponds to the i-th diagram in fig. 2, and 

Of" = \n c C f (2C f - \c A ) = ~ Cf" = \n c C f 1 -C a = 3. (10) 

The second term in eq. (]?]) is characteristic of non-Abelian theories and would be absent in 
any Abelian model. 

Finally, for the production of two quarks, a gluon and a photon (|3|), the matrix element 
squared is (up to a constant factor, see later on): 

\M\ 2 = C^\e q M a + e e M ISR \ 2 , (11) 

with 

C^m = N C C F , (12) 

where M a is the sum of the diagrams of fig. 2a, Misr the sum of the diagrams of fig. 2c, and 
e e and e q are the electric charges of the electron e and of the quark q in the final state. 

The Abelian model we compare with QCD is the one introduced in [26] J Here gluons 



have no colour and no self-coupling: therefore, only diagrams of fig. 1 and fig. 2a survive. The 
cross sections of processes (|I])-(^D for the Abelian case can be obtained from the QCD ones by 
simply replacing the group constants of QCD by those appropriate for the Abelian model: i.e., 
Ca = 3— >0,CV = 4/3^1 and T F = 1/2 — > 3. The 'Abelian coupling constant" a a is fixed 
to (4/3)a s , so that the ratio of the two-jet and three-jet cross sections agrees with experiment. 

We have analyzed the processes (|T])-([3]) adopting four different jet-finding algorithms. They 
are identified through their clustering variable y^. They are the JADE scheme (J) |28| based 
on the variable 

2EiEAl - cos 9, 



= -^v- ^ (13) 



and its "E" variation (E)Q 
the Durham scheme (D) B 



.E (Pi+Pj) ■ (Pi+Pj) 



V* = ^ * J ' ^ , (14) 



d 2mm(E 2 ,E 2 )(l-cos6 lJ ) 
Vii = , ( 15 ) 



4 Of course, we know that such a model has been already ruled out in other contexts, e.g., by measurements 
of the energy dependence of multi-jet production rates in e + e~ annihilation f27j], but we regard it mainly as a 
useful tool to demonstrate the sensitivity of the introduced angular variables to the various features of QCD. 

5 At lowest order, the E and JADE schemes are equivalent for massless particles. 
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and the Geneva algorithm (G) J3(J 



G = S EjEjjl -COS%; 

• A ' ; 9 (/•;, + /•;/•' ' 



For all of them the two (pseudo)particles i and j (with energy Ei and Ej, respectively) for 
which i/ij is minimum are combined into a single pseudoparticle k of momentum given by 
the formula 

P k = P i + P j . (17) 

The procedure is iterated until all pseudoparticle pairs satisfy y^ > y cut . The various charac- 
teristics of these algorithms are well summarized in ref. [[I(J. In our lowest order calculation, 
the four jet cross section for a given algorithm is simply equal to the four parton cross section 
with a cut yij > y cut on all pairs of partons 

Concerning the numerical part of our work, we have taken a em = 1/128 and sin 2 9w = 
s^r = 0.23, while for the Z° boson mass and width we have adopted the values M^o = 91.1 
GeV and T^o = 2.5 GeV, respectively. For the quarks we have: m c = 1.7 GeV and m;, = 5.0 
GeV while the flavours u, d and s have been considered massless. Finally, the strong coupling 
constant ot s has been set equal to 0.115. 



3 Angular variables 

We study the following four variables: the modified Nachtmann-Reiter angle, 9* NRl the Beng- 
tsson-Zerwas angle, Xbz, a modification of the Korner-Schierholz-Willrodt angle we denote by 
&*kswi an d the angle between jets 3 and 4, # 34 , in two different situations: (a) with 6-tagging 
and (b) without 5-tagging. 

In case (a), we consider four-jet events where two of the jets contain a b or a b. Let us 
call them jet 1 and jet 2. It does not matter which one is the b- and which the 6-jet. We do 
not make any assumptions about jets 3 and 4; they may be either gluon or quark jets, or even 
6-jets. Then, in terms of the three-momenta p±, . . . ,p4 of jets 1, ... ,4, the angles 6* NRl xbz 
and #34 are defined by 

e NR = l (Pi ~P2,P3 -Pi), (18) 
Xbz = L{p x x p 2 ,p 3 x p 4 ), (19) 

and 

9 u = l(p 3 ,p 4 ). (20) 

For events where 

Wi +P3\ > \pi +P±\ (21) 

we define 

®ksw = x P^P2 x n)- (22) 
In the opposite case, we define &* KS w with p 3 and p 4 interchanged. The definition in eqs. (pl"l)- 



(22) is equivalent to the original definition of $^5vk[13| in events where the thrust axis is along 



Pi +p 3 or p x +p 4 . 
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In situation (b), where there is no 6-tagging, we label the jets according to their energy, 
such that Ei > E 2 > E 3 > E4, and then define the angles by eqs. (|l"8l)-([2"2j) as before. 

By considering the polarization of the gluon in e + e _ — > qqg and the final state helicities in 
the subsequent splitting g — > gg or g — > q'q', one finds that the e + e~ — ► gggg cross section is 
concentrated near cos9^ R « ±1, whereas in e + e~ — > qqq'q', the cross section is largest around 
cos6nr ~ 0. In the case of the Bengtsson-Zerwas angle, one expects the g — > gg contribution 
to be rather flat in the corresponding distribution if compared with the g — > qq one, which 
generally peaks at xbz ~ 90°. The original Korner-Schierholz-Willrodt angle is defined for 
events for which there are two jets in both the hemispheres separated by the plane perpendicular 
to the thrust axis: it is the angle between the oriented normals to the plane containing the 
jets in one hemisphere and to the plane defined by the two other jets. The advantage of the 
modified definition &* KSW adopted here is that it allows us to include the complete 4jet-sample 
in the analysis, without having to discard events with three vectors in the same hemisphere. 
In the splitting process g — > gg, the two planes tend to be parallel, with the two gluons on 
the same side, i.e., Q^sw ~ 71 1 whereas the planes are preferentially orthogonal for g — > qq. 
Finally, gluons from the triple gluon vertex g —>■ gg and the second pair of quarks from g — > qq 
are expected to be closer together than gluons from double bremsstrahlung, and this should be 
evident by looking at the angle between the two softest jets. More details on these arguments 



are given in ref. 14 



The results for the Nachtmann-Reiter angle are shown in figs. 3a & b. In the case of 
6-tagging (fig. 3a), the distributions are even functions of cos9^ R , because replacing 6* NR by 
71 — 9* NR is equivalent to interchanging p 3 and pk, and we do not make any distinction between 
jets 3 and 4. The fefegg-distributions have peaks at cosB^r = ±1, whereas the maxima of the 
bbqq are at cos8^ R = 0. We also note that the fefegg-distributions in QCD and in the Abelian 
model are different. Without 6-tagging, the distributions of cos 9^ R (fig. 3b) are skewed to the 
left. The reason for this asymmetry is kinematical and can be understood by considering events 
where all four jets are all close to one common axis (with small angles between them in order to 
satisfy the y-cuts). Then, by momentum conservation, the most and the next most energetic 
jets, 1 and 2, must go in opposite directions. If we further restrict our attention to events where 
jets 3 and 4 also go in opposite directions, then energy-ordering and momentum conservation 
together imply that jet 3 must be parallel to jet 2, which gives 6 NR = ir, and cannot be 
parallel to jet 1, which would give 6 NR = 0. Although the qqgg and ggg'g'-distributions are 
still different, the difference between QCD and the Abelian model in the gggg-distribution is 
washed away. 

In figs. 4a & b, the distributions are shown of the Bengtsson-Zerwas angle, with and without 
6-tagging, respectively. With 6-tagging the distributions are again symmetric, for the same 
reason as in the case of Nachtmann-Reiter. At xbz = tt/2, there is a peak in the bbqq- 
distributions and a dip in the 66gg-distributions. Without 6-tagging, the fefegg-distributions 
are shifted to lower values of Xbz, whereas the fefegg-distributions are shifted slightly towards 
higher \bz- 

The distributions where the rewards for 6-tagging are largest are probably those of the angle 
®*kswi shown in figs. 5a & b. With 6-tagging, the differences between the bbgg distributions 
in QCD and the Abelian model become even more clear than in xbz- In this respect the 
modified definition of &*kswi initially adopted in order to avoid loss of statistics, turns out to 
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be extremely successful. 

Finally, figs. 6a & b show the distribution of cos# 34 . Here we see, as expected, a tendency 
for the gluons from the triple gluon vertex, and the quarks qq in hhqq events, to be closer 
together than the gluons in the Abelian model. In fig. 6a, there is also a large concentration of 
bbqq events near cos #34 — — 1. They come from the region of phase space where the 6-quarks 
are relatively soft, where the cross section is dominated by the diagrams in which the Z° is 
coupled directly to the quarks qq. That explains why the peak disappears completely when 
the jets are energy-ordered (fig. 6b). 

We should warn the reader that the value of y cu t we have chosen for using with the Geneva 
algorithm corresponds to a looser cut than the ones we use with the other algorithms. This 
is because, unlike the other jet defining variables, the definition ( |i~6"D of y G does not contain s 
explicitly, and therefore, the Geneva algorithm allows the energies of the partons to be much 
smaller for a given value of y cut . As a result, we can get very close to the singularities of the 
matrix elements, where we expect radiative corrections to be large (figs. 5 and 6). The large 
peaks in the cross section also make it more difficult to integrate by VEGAS, as can be seen 
from the statistical fluctuations in the distributions of figs. 3 & 4. Therefore, a larger value of 



Vent wou ld be needed to obtain reliable predictions from our tree- level calculation |30| , but we 
prefer to keep the value shown just to illustrate what happens. 



4 Photon sample 

In the literature [|T^, it has been argued that one might be able to see the non- Abelian 
structure of QCD by comparing the distributions of four-jet events with those of events where 
three jets and a photon are produced. Since the diagrams of fig. 2a are the same, up to a 
constant factor, for the processes (§) and @, any differences in their distributions must be 
due to the non-Abelian diagrams of fig. 2b and the ISR diagrams of fig. 2c. As we shall show 
below, the contribution of the latter can be made negligibly small by applying suitable cuts. 
If one then assumes that the four-jet events are predominantly qqgg events, as is true in QCD, 
differences between the four-jet and the three-jet plus photon distributions can be regarded as 
evidence for the non-Abelian contribution |.Mb| 2 in eq. (^). 

What effect would 6-tagging have on such an analysis? Presumably, the distinctions be- 
tween the distributions would become more clear, but the number of events would be smaller. 
Moreover, selecting events with a 66-pair in the final state increases the relative number of 
unwanted four-quark events, since five of the fifteen flavour combinations qqq'q', but only one 
of the five combinations qqgg (where q, q' = u,d,c,s,b), contain at least one 66-pair. It is 
possible to reduce this contamination of the four-jet sample by imposing cuts, at the cost of 
a further loss of statistics. The crucial question is, whether, in the end, the event rates would 
still be large enough to allow a study of the distributions. 

The ISR diagrams are important for photons that are either soft or nearly collinear with the 
incoming electrons and positrons. Therefore, we can eliminate them by imposing cuts on the 
photon energy and on the angle Qheam-i between the photon trajectory and the e 1 * 1 beams. 
Since we wish to compare the sample of three-jet plus photon events with an equivalent sample 
of four-jet events, we must treat both samples on exactly the same footing. In the four-jet 
sample we impose the same cuts on all four jets, or at least on the two non-6 jets, in case we 
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select events where two of the jets are tagged as 6-jets. This implies that in the three-jet plus 
photon sample, we must also impose the same cuts on the jet energies Ej Gt and their angles 
with respect to the beams 9b ea m-jet- In the 6-tagged case, we only apply these cuts to the non-6 
jet. 

It turns out that demanding that |cos#6 e am-7l < 0-9 an d -E7 > 10 GeV is sufficient to 
remove the effect of the ISR diagrams. This is illustrated, for example, in fig. 7, where we 
have implemented the cut | cos Obeam-j 1 < 0.9 and plotted the differential distribution in the 
photon energy of the three-jet plus photon cross section twice: taking the ISR diagrams into 
account and omitting them. Above our 10 GeV cut on E 1 , the two distributions are the same. 
We have checked that, with these cuts, the distributions of other variables also look the same 
with and without ISR. 

The total cross section for process (|^) is given in tab. I for several values of y cut . The effect 
of the ISR is never greater than a few percent. 

In tab. II we show the qqgg and qqq'q' components of the total e + e~ — > 4jet cross section. In 
tab. Ha, where we assume two 6-jets are tagged, we have imposed two additional cuts, namely 
I cos 6* NR \ > 0.5 and xbz < 50° or xbz > 130°. As can be seen in figs. 3 & 4, this reduces the 
relative size of the qqq'q' component, making it less than about 8% of the total four jet cross 
section for all values of y cut shown in the table. We have not imposed these cuts in tab. lib, 
which shows the results if one does not select 6-tagged jets, because there the qqq'q' component 
is already quite small without them. 

A comparison between the e + e~ — > qqgg and e + e~ — > qqg'y cross sections is made in fig. 8. 
They are displayed as a function of y cut for each of the jet-finding algorithms. The dotted curves 
marked "real" show the actual qqg^i cross section. The dashed curves marked "renormalised" 
show the qqg'j cross section multiplied by 

(23) 



Apart from the small ISR effect, this is exactly the contribution of the QED-like graphs, i.e., the 
term |7W a | 2 in eq. ([?]), to the e + e~ — > qqgg cross section. The factor of two in the denominator 
is the symmetry factor needed to account for the two identical gluons in the qqgg final state. 
The non-Abelian contribution |.Mb| 2 is the difference between this "renormalised" qqg'j cross 
section and the qqgg cross section. 



5 Mass effects 

In this section, we examine the numerical importance of taking the 6-quark mass into account 
exactly, as we have done in all the calculations we have discussed until now. The 6-quark 
mass does not only affect the total cross section |22], |23], [24|, but also some of the angular 



distributions. In figs. 9, 10 & 11, we present plots illustrating this for the processes e + e~ — > 
bbgg, e + e~ — > bbuu, and e + e" — ► bbbb. We do not show, e.g., e + e~ — > bbdd, because the results 
are quite similar to those of e + e~ — > bbuu. In each case, we compare the distributions obtained 
for m& = 5 GeV with the ones we would find if we neglected m&. 

The differences turn out to be fairly small in the bbg g-process, but in the other ones they 
are quite large, particularly in the distributions of xbz and &* KS w This is true both for the 
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distributions where the fo-jets are identified, see figs. 10a & 11a, and for the distributions where 
the jets are simply ordered in energy, figs. 10b and lib. These mass effects also depend on the 
particular jet-finding algorithm used. As an example we plotted the case of the E scheme, for 
which the mass effects are larger. 

In the equal flavour process, e + e _ — > bbbb, there are three different combinations of b- 
(anti)quarks that can be tagged: two 6's, one b and one b, or two 6's. That is why fig. 11a 
shows two sets of curves, one for the bb case and one for the bb casef]. We also note that even in 
the bb case, the distribution of cos 9* NR is very different from the corresponding distribution in 
the unequal flavour process. Due to the gluon propagators in the last four diagrams of fig. 1, it 
peaks at co$6* NR = ±1, making it look similar to the cos distribution of the bbgg process. 



6 Conclusions 

We have given our results in terms of total and differential cross sections. To convert them 
into event rates, they must be multiplied by the luminosity and the efficiency for tagging two 
6-quarks. We have not done this because the numbers differ from one experiment to another, 
and could still change as the techniques are improved. However, to get an idea of what one 
might expect, let us suppose the integrated luminosity is 100 pb _1 and the tagging efficiency 
50%. Then there will be roughly 4 x 10 4 four jet events with two tagged fe-quarks (using the 
JADE scheme with y cut = 0.01, as in figs. 3-6). This is an order of magnitude less than the 
total number of four jet events, but in return we gain greater power to discriminate between 
the various terms in the four jet cross section, particularly between the Abelian, QED-like 
term and the non-Abelian term in the two quark, two gluon cross section. 

In a realistic analysis, one should take the probability of misidentifying other particles as 
6-quarks into account. We do not expect c-jets mistagged as 6-jets to distort the distributions 
severely, although they would increase the total number of events, because they would still 
be correctly classified as quark jets. It is, of course, important that the number of gluon jets 
tagged as 6's be as small as possible. The probability of a 6-quark being created during the 
fragmentation of a gluon or a lighter quark is believed to be negligible W§ . 



The number of three jet 7 events is reduced even more severely by demanding two tagged 
6-quarks, because of their small electric charge (—1/3 instead of 2/3). With the luminosity 
and efficiency given above, and the cuts of fig. 8a, one would expect of the order of 50 events, 
which is not enough to study any distributions. So, in the photon sample, it is better to keep 
all events. 

Finally, we indicated the effects of the 6-quark mass. Normally, they can safely be neglected 
because only a small fraction of the four jet events contain 6-quarks, and moreover, most of 
those are bbgg events, where the effects are small. This is no longer true with 6-tagging, because 
now all events contain at least two massive 6-quarks, and the relative number of four quark 
events is higher, and the mass effects are especially important when comparing QCD with the 
Abelian model, where as many as 30% of the events are four quark events. 



6 In figs. 3a, 4a, 5a & 6a, we assumed a b and a b were tagged. If that assumption is not true, the contribution 
of the e + e~ — > bbbb events to those distributions has to be replaced with a weighted average of the two sets 
of curves in fig. 11a. However, since the contribution of the e + e~ — > bbbb events is small, the distributions in 
figs. 3a, 4a, 5a & 6a would hardly change. 
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Table Captions 

tab. I Cross sections in picobarns of the processes (a) e + e~ — ► bbg^y (i.e., 6-tagging) and (b) 
e + e~ — > J2 q <l<igj (i-e., no 6-tagging), with and without ISR, for three different values of 
y cut of each jet-finding algorithm, using the notation (a IS R, cr no iSR\ ycut)- 

The following additional cuts have been implemented: | cos #beam-7, 9 1 < 0.9 and E 7>g > 10 
GeV, for case (a), and | cos^ b eam-7, s ,, 3 ,g| < ^-9 and -E 7 , s , g ,g > 10 GeV, for case (b). 

tab. II Cross sections in picobarns of the processes (a) e + e~ — > bbgg (x) and e + e" — > J2 q bbqq 
(y) (i.e., 6-tagging) and (b) e+e~ -> Y, q W99 (x) and e+e~ ->■ Y. q Y. q >qw'q' (y) (i-e., no 
6-tagging), for three different values of y cut (z) of each jet-finding algorithm, adopting the 
notation (x, y; z). 

The following additional cuts have been implemented: | cos0beam-g| < 0.9, E g > 10 GeV, 
\cos9^ R \ > 0.5 and xbz < 50° or xbz > 130°, for case (a), and | cos9h ea ,m- g , q , q \ < 0.9 
and E g ^ q £ > 10 GeV, for case (b). 



Figure Captions 

fig. 1 Feynman diagrams contributing in lowest order to e + e~ — > qqq'q', where q^' = u, d, s, c 
and 6. If q 7^ g' only the first four diagrams contribute. The internal wavy line represents 
a photon or a Z°. The particles are labelled as in eq. ([[]). 

fig. 2 Feynman diagrams contributing in lowest order to e + e~ — > qqgg (a and b) and e + e~ — ► 
gg<?7 (a and c), where q = u,d, s,c and 6. The internal wavy line represents a photon or 
a Z°, while the external jagged line represents a gluon or a photon, as appropriate. The 
particles are labelled as in eqs. (0)-(0). 

fig. 3 Distributions in the cosine of the modified Nachtmann-Reiter angle, cos$yy^>, in (a) 
e+e~ -> bbgg,J2 q bbqq~ (i.e., 6-tagging) and (b) e+e~ -> E q qqgg,E q E q > qqq'q' (i.e., no 
6-tagging), for the various jet-finding algorithms. 

fig. 4 Distributions in the Bengtsson-Zerwas angle, Xbzj m ( a ) e + e _ — > bbgg,J2 q bbqq (i.e., 
6-tagging) and (b) e + e~ — > J2 q qqgg :J2 q J2 q ' qqq'q' (i- e - ; 110 ^-tagging), f° r the various 
jet-finding algorithms. 

fig. 5 Distributions in the modified Korner-Schierholz-Willrodt angle, &*ksw-i m ( a ) e+e ~ ~> 
bbgg,E q bbqq (i.e., 6-tagging) and (b) e+e~ -> E q qqgg,E q E q > qqq'q' (i.e., no 6-tagging), 
for the various jet-finding algorithms. 

fig. 6 Distributions in the cosine of the angle between the vectors ^3 and P4, cos $34, in (a) 
e+e" -> bbgg,J2 q bbqq (i.e., 6-tagging) and (b) e+e~ -> E q qqgg,E q E q > qqq'q' (i.e., no 
6-tagging), for the various jet-finding algorithms. 

fig. 7 Distributions in energy of the photon in (a) e + e~ — > 66(77 (i-e., 6-tagging), with 
I cos6 l beam-7,g| < 0.9 and E 7t9 > 1.0 GeV, and (b) e + e~ — > J2 q qqgi (i-e., no 6-tagging), 
with | cos #beam-7,g,g,g| < 0.9 and ^ 7 ,g, g ,f > 1-0 GeV, with and without ISR, for the various 
jet-finding algorithms. 

fig. 8 Cross sections of the processes (a) e + e~ — ► 66(7(7, 66(77 (the latter both real and renor- 
malized to the former, see in the text) (i.e., 6-tagging), with | cos 6 l beam _ 7j9 1 < 0.9 and 
E J>g > 10 GeV, and (b) e + e~ — > J2 q qqgg,12 q qqgi (the latter both real and renormal- 
ized to the former, see in the text) (i.e., no 6-tagging), with | cos9beaia--y,g,q,q\ < 0.9 and 
E 1: g, q ,q > 10 GeV, as a function of y cut for the various jet-finding algorithms. The follow- 
ing angular cuts have been also implemented in case (a): | cos > 0.5 and xbz < 50° 
or X BZ>130°. 

fig. 9 Mass effects in the angular distributions of the process e + e~ — > bbgg using the E jet 
finding algorithm. In (a) 6-tagging is assumed, in (b) the jets are energy-ordered. The 
curves denoted by bbgg are for massive, and those denoted by ddgg for massless 6-quarks. 

fig. 10 Mass effects in the angular distributions of the process e + e~ — > bbuu (q 7^ 6), using the 
E jet-finding algorithm. In (a) 6-tagging is assumed, in (b) the jets are energy-ordered. 
The curves denoted by bbqq are for massive, and those denoted by ddqq for massless 
6-quarks. 



11 Mass effects in the angular distributions of the process e + e~ — > bbbb, using the E jet- 
finding algorithm, (a), solid and short-dashed lines: one b and one b are tagged; dotted 
and long-dashed lines: two fe's are tagged. In (b), the jets are energy-ordered. The curves 
denoted by bbbb are for massive, and those denoted by dddd for massless 6-quarks. 
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